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Transport of molecular motors, stimulated by interactions with specific links be- 
tween consecutive binding sites (called "bridges" ) , is investigated theoretically by an- 
alyzing discrete-state stochastic "burnt-bridge" models. When an unbiased diffusing 
particle crosses the bridge, the link can be destroyed ("burned") with a probability 
p, creating a biased directed motion for the particle. It is shown that for probability 
of burning p = 1 the system can be mapped into one-dimensional single-particle hop- 
ping model along the periodic infinite lattice that allows one to calculate exactly all 
dynamic properties. For general case of p < la new theoretical method is developed, 
and dynamic properties are computed explicitly. Discrete-time and continuous-time 
dynamics, periodic and random distribution of bridges and different burning dynam- 
ics are analyzed and compared. Theoretical predictions are supported by extensive 
Monte Carlo computer simulations. Theoretical results are applied for analysis of 
the experiments on collagenase motor proteins. 

INTRODUCTION 



ar motors, or motor proteins, play important role in the functioning of biological 
. They transform a chemical energy into mechanical work with a high ef- 
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ficiency and move along linear molecular tracks. Conventional molecular motors fuel their 
motion by hydrolyzing adenosine triphosphate (ATP) or related compounds, and the chem- 
ical energy influx creates a biased directed motion along one-dimensional protein filaments. 
During the translocation of motor proteins the chemical state of linear protein tracks remains 
unchanged. 

Recently, a new type of molecular motor, independent of ATP hydrolysis, has been dis- 



covered jj,|5|. A special protein called collagenase moves along collagen fibrils, which are 



the main component of extracellular matrix 



l|, and it cleaves the filament at special sites. 



Collagen proteolysis provides a chemical energy that produces a directed motion of the mo- 
tor protein molecule. It was found that the random motion of enzyme molecule can be 
transformed into the biased diffusion along the linear track if, after the cleavage of collagen, 
the motor protein could only move on one side of the broken link jsj]. 

It was suggested that the dynamics of collagenases can be explained by a so called "burnt- 
bridge" model (BBM) DDDQ . According to this approach, the motor protein is viewed 
as a random walker hopping along the one-dimensional lattice that consists of two types 
of links: strong and weak. Strong links are not affected when the particle passes them, 
however crossing the weak links might destroy them with a probability < p < 1, and 
the random walker is not allowed to move again over the "burnt" links. Using a spatial 
continuum approximation, that neglects the underlying lattice, the burnt-bridge model has 
been discussed in the limiting cases of very low p and p — > 1 [6j. This approximation is 
only valid in the limit of very low concentrations of the special links, although in biological 
systems the density of bridges might be significant jj, 0]. Antal and Krapivsky [2] have 
investigated BBM by using a discrete-time random walk approach that allowed to calculate 
exactly velocities and dispersions for some sets of parameters for several versions of the 
model. However, discrete-time dynamics is unrealistic for motor proteins since their dynam- 
ics is tightly coupled to a set of chemical transitions -8], and time intervals between these 
transitions are distributed exponentially according to Poisson statistics. A discrete-time 
statistics implicitly assumes a delta-function waiting time distribution between consecutive 
steps. In addition, the approach presented in Ref. 7J is mathematically very involved, and 
not all relevant dynamic properties have been computed. The goal of this paper is to in- 
vestigate more realistic continuous-time discrete-state stochastic burnt-bridge models using 
simple analytical approaches and extensive Monte Carlo computer simulations. 



MODEL 



Let us consider a single random walker that moves along an infinite one-dimensional 
lattice as shown in Fig. 1. The particle can hop in both directions with the same rate 
vq = 1, i.e., initially there is no bias in the motion of the random walker. We set the size 



3 



of the lattice spacing to be equal to one. There are two types of links between consecutive 
binding sites on the lattice. The random walker does not interact with the strong links, while 
crossing the weak link leads to burning of the bridge with the probability p. We assume 
that after the burning of the bridge the particle is always on the right side from it, and there 
are no burnt bridges at the initial time. Thus the particle is never trapped between broken 
links, and it moves continuously forward. 

We distinguish two different possibilities of burning the bridge. In the first option, the 
bridge is burned only when the particle crosses it from left to right, but nothing happens 
when the walker moves from right to left. This model is called a forward BBM. In the 
second option, the particle burns the bridge on any forward and backward motion across 
the weak link, describing a forward-backward BBM. It is important to note that both cases 
are identical for p — 1, while for p < 1 dynamics differs significantly for different burning 
scenarios. Dynamic properties of the system also depend on the details the particle behavior 
near the already burned bridge []|. In the standard version of BBM the particle sitting next 
the broken link can only move forward to the right direction. In the modified version of 
BBM the walker tries to cross the burnt bridge, but the attempt fails and the particle does 
not change its position. 

The dynamics of the motor protein molecule in different BBMs is directly determined by 
two parameters: the probability of burning p and the concentration of weak links c. The 
distribution of bridges strongly influences the dynamic properties \\ , and we will consider 
two different cases: periodic and random distributions. The case when the bridges are found 
at equal distance N — 1/c lattice spacings apart from each other describes the periodic 
distribution, while in the random case we postulate that each link is a bridge with the 
probability c every time the walker visits it. Note, however, that other realizations of the 
random distribution of bridges are also possible. 



BURNT-BRIDGE MODELS WITH p = 1 

First, we consider a special case of BBM with p = 1 when any crossing of the bridge leads 
to its burning. In addition, it is assumed that the first burnt bridge was crossed from the 
left. Thus the particle drifts to the right (forward) direction. For the periodic distribution 
of bridges the dynamic properties of the system can be evaluated explicitly by mapping the 
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model into the motion of the random walker on infinite one-dimensional periodic lattice, for 
which exact sciaticas have been derive, by Derrida to 198 3 flflQ. The particle at site 3 
(j — 0, 1, • • • , N — 1) moves forward (backward) with the rate Uj (wj), and the periodicity 
requires that Uj = Uj+N and Wj = Wj + N, where N is the size of the period. For this model 
the explicit expression for the mean velocity V is given by 
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while the dispersion, or diffusion coefficient, can be found from 
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The BBM with p = 1 corresponds to the particle hopping model with 



(3) 



(4) 



1 for all 



j with the exception of Wq = 0, that reflects the complete burning of the bridge. Then from 
Eq. (JIJ the explicit formula for the velocity is given by 

2 2c 



V 



(5) 



iV + 1 c+1 
The dispersion can be obtained from Eq. yielding 

3(c+l)2 w 

The dynamic properties of BBM with p = 1 are shown in Fig. 2 as functions of the 
concentration of bridges. The velocity is a monotonically growing function of c, while the 
dispersion is always decreasing. It is interesting to consider their behavior in limiting cases. 
When every link on the lattice is a bridge (c = 1) we obtain V — 1 and D — 1/2, as expected. 
More interesting is the limit of c — > 0. The velocity goes to zero, V ~ 2c, as expected, and 
this reproduces the result of the continuum theory of Mai et al. (|. But, surprisingly, it 
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is found that D BBM (c = 0) = 2/3 and it is not equal to the dispersion of a free particle 
on the lattice without bridges, Df ree = 1. This jump in the dispersion might correspond 
to a dynamic transition, indicating qualitatively different behavior of the particle moving 
in the system without weak links andparticle diffusing in BBM. Similar behavior has been 
observed in the discrete-time BBM 7J. However, this interesting phenomenon requires a 
more careful and detailed investigation. 

The mapping of BBM to the particle hopping model on the periodic lattice can also be 
used to investigate discrete-time burnt-bridge models. It was found [l| that in this case 
the explicit expression for the velocity is the same as in the continuous case, while for the 
diffusion coefficient it can be shown that Ddiscr = D cont — V 2 /2, and D cont and V can 
be obtained from Eqs. and (J2J). However, the rate constants Uj and Wj have to be 
reinterpreted as jumping probabilities, leading to 

Uj = Wj = ^ forj = l,2,--- ,N-1. (7) 

At the same time, for the particle near the burned bridge the rates are uq = 1 and wq = 
for the standard BBM, while Uq = 1/2 and Wq = for the modified BBM. Then the dynamic 
properties of the standard discrete-time BBM are 

V = c, D= 1 -(l-c 2 ), (8) 

and for the modified BBM we obtain 

c 1 c 2 

V = 7+l , D= 3(c+l) ~ 6(c+l) 2 ' (9) 

These results fully agree with those obtained in Ref. [tJ for periodic bridge distribution. 
Note, however, that our approach based on Derrida's analysis is much more straightforward. 



BURNT-BRIDGE MODELS WITH p < 1 

Derrida's method cannot be easily extended to include the general case of bridge burning 
with the probability p < 1 since it is not clear how the hopping rates Uj and Wj depend 
on p. To calculate dynamic properties of BBM we introduce a new method that is very 
general and applicable for many different systems. Below we analyze in detail several BBM 
to illustrate our approach. 
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Continuous-time forward BBM for the periodic distribution of bridges 

Now consider a continuous-time BBM with the particle destroying the link with the 
probability p only when the bridge is crossed from left to right. The distribution of weak 
links is periodic, i.e., they are positioned N — 1/c sites apart from each other. Let us define 
Rj{t) as the probability to find the random walker j sites apart from the last burnt bridge at 
time t. The probabilities Rj(t) are associated with the moving system of coordinates with 
the last burned bridge always at the origin. The corresponding kinetic scheme is shown in 
Fig. 3a. The dynamics of the system can be described by a set of Master equations. For 
the lattice sites not connected to weak links or connected only from the left side we have 

kN ^ 1 ^ ' = RkN+i-l{t) + RkN+i+l{t) — 2RkN+i{t), (10) 

for k — 0, 1, 2, - • • and i — 1, 2, • • • , N — 1. For all other lattice sites (except the origin) one 
can show that 

dR ™® = (1 - p)R kN . l (t) + R kN+ i(t) - 2R kN (t) for k = 1, 2, • • • . (11) 
The dynamics is different at the origin (see Figs. 1 and 3a) where 

= pf2[R kN -i(t)) + Ri(t) - R (t). (12) 

k=l 

In addition, there is a normalization condition, 

oo N-l 

5^iWW = l, (13) 

k=0 i=0 

which is valid at all times. 

In the stationary-state limit (t — > oo) we have f t Rj(t) = 0, and Eqs. (HUJ), (HQ) and flUJ) 
simplify significantly, 

oo 

Ri = R -ipS S = Y^ RkN-i, for i = 0,1, ■ ■ ■ ,N; (14) 

k=l 

and 

RkN+i = (i + l)RkN - «(1 - p)RkN-i for z = 0, 1, • • ■ ,N; k = 1, 2, • • • . (15) 
To find the solutions of Eqs. ()14j) and (|15jl we assume that it has the following form 

R kN+l = R e ak -tB(k), (16) 
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where a and B(k) are unknown parameters. This expression assumes a linear dependence 
between the weak links and an exponential decrease for consecutive periods. To eliminate 
B(k), we observe that Eqs. (fUjl and (fTHjl simultaneously hold for i = N, giving us 

R( k +i)N = i? e a(/c+1) = R k N + N = Roe ak - NB(k), (17) 

that leads to 

B{k) = ^e ak {l-e a ). (18) 
To obtain the expression for a, we substitute RkN+i from Eq. ()16|) in the recurrent formula 

m, 

R e ak - iB(k) = (i + l)i2 e a * - z(l - p) [Roe ^ - (JV - l)fl(Jfe - 1)] . (19) 

Then using the expression (fTHj) for -B(/c) and defining x = e a , Eq. (|T9"j) eventually reduces 
to 

x 2 - [2 + p(N- l))x+ (1 -p) = 0. (20) 

Since RkN+i is a decreasing function, the parameter a is always less than one, and we choose 
a physically reasonable solution of Eq. (|2U|) , 

x = 1 + ip(iV - 1) - ^Vp 2 ( n ~ !) 2 + 4p/V- (21) 

Note that generally < x < 1 and x = for p = 1, while x = 1 for p = 0. The probability to 
find the particle next to the burnt bridge, Rq, can be found from the normalization condition 
(JT3J) and Eq. J2H, yielding 

R = 2(1 -x) = 2p = 2pc 

N + l + x(N-l) y/p*(N - l) 2 + 4pN v/p 2 (l-c) 2 + 4pc' 

and the solution for RkN+i in terms of x and Rq is finally given as 



RkN+i — RqX 1 



1 -£(!-*> 



(23) 



Exact solutions for the probabilities to find the particle at given distance from the 
last burnt bridge allow us to calculate many dynamic properties of the system. Specifically, 
the mean velocity of the particle is given by 

oo 

V = J2(u j -w j )R j . (24) 

j=0 
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However, since the rates to go forward and backward are the same at all sites except the 
origin, all contributions to the velocity cancel out and only j — term survives, 



For p = 1 this equation reduces to V(c,p = 1) = 2c/ (c + 1), which agrees exactly with the 
result © derived using the mapping to Derrida's random hopping model. In the limit of 
p — > and c (p < c) Eq. (|23j) becomes V ~ ^Jpc as compared to V ~ \J1pc in 
The discrepancy is due to the fact that in our model only forward burning is considered, 
while in Ref. (p) the bridge can be burned on both forward and backward motions of the 
motor protein molecule. The expression for the velocity (J25)) also simplifies in the case of 
N = c = 1, giving V = ^Jp. 

Another useful dynamic property that can be obtained from our analysis is average and 
average-squared distances between the particle and the last burnt bridge, namely, 

oo oo 

(j) = J2jRv and (f) = J2fR r (26) 

j=0 j=0 

Substituting Eq. (f2*??j) into the expressions (f2T)j) we obtain 

... „ f N 2 x N 2 -l) 



,N 3 x(x + l) 1 2 + l 



and 

where x and i?o are given by Eqs. (|21|) and (|22j) . For p = 1 one can show that 

(i) = (N - l)/3, (j 2 ) = iV(iV-l)/6. (29) 
In the limit of p C c we have 

(j) ~ v^V>, (j 2 ) ~ 2iV/p, (30) 
while in the opposite limit of c <C p the results are 

(j) ~ iV/3, (j 2 ) ~ iV 2 /6. (31) 
For the special case of iV = c = 1 we obtain 
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We can also calculate the average and average-squared distances between the consecutive 
burnt bridges. Because the burning is a stochastic process, i.e., p < 1, this distance generally 
deviates from N. In order to compute these properties we introduce PkN as a probability 
that the next bridge that burned is kN sites from the given burnt bridge. Then 

^ pRkN-i 



kN 



A; = 1, 2, - - - . (33) 



The average and average-squared distances between the consecutive burnt bridges can be 
found from the following expressions 

oo oo 

(/) = kNP kN , (I 2 ) = J2( kN ) 2p kN- (34) 

k=l k=l 

By substituting Eqs. (J23|) and (J3H|) into and using Eq. (J2T)) we found that 

W " <««> = ^f>. (35) 
1 — x (1 — xj 2 

The average distance (/) is related to the fraction of unburnt bridges q via the following 
expression 

q= mg_ii =1 . m)=x , (36) 

which gives the physical meaning of the variable x used in our analysis. In the limiting cases 
we obtain the following results: 

(/) = N, (I 2 ) = N 2 , for p = 1; 

(l) = l/y/p, (l 2 ) = 2/p-l/^, for N = c= I; 



(I) = V /]V>, (I 2 ) = 2N/p, for p < c; 

(/) = jV, (/ 2 ) = A^ 2 , for p > c. 

Discrete-time forward-backward BBM for the periodic distribution of bridges 

Our theoretical method is very general and it can be applied to any version of BBM. 
As another example, we analyze a motion of the particle in discrete-time framework with 
periodically located bridges that can burn with the probability p during the forward and 
backward transitions. We again introduce a modified kinetic scheme by putting the last 
burnt bridge at the origin, as shown in Fig 3b. The important difference between continuous- 
time and discrete-time models is the fact that instead of transition rates in the continuous- 
time case transition probabilities are used in the discrete-time dynamics. As a result, the sum 
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of transition probabilities at each site is equal to unity. We assume that at the sites that do 
not support the weak links the transition probabilities are equal, i.e., u k N+i = w k N+i = 1/2 
and WkN~i = UkN = 1/2 for any k and for i — 1, 2, • • • , N — 2. The transition probabilities 
are also equal at the sites that surround the weak links, giving uun-i = WkN = (1 — p)/2 for 
k — 1, 2, • - • (see Fig. 3b). In addition, at the origin it is assumed that uq = 1 and wq = 0. 

As before, we define Rj(t) as the probability to find the random walker at the position j 
at time t. The dynamics of the system is governed by a set of discrete-time Master equations 
that are different from the continuous-time version [?], E|. For the lattice sites that do not 
support the weak links the corresponding Master equations are 

RkN+i(t + 1) = ~RkN+i-l(t) + -RkN+i+l(t) > (38) 

for any k and with i — 1,2, • ■ ■ ,N — 2 (with the exception of Rq and R\ — see below). For 
the lattice sites around the bridges we have 

1 1 -p 

RkN-l(t + 1) = -RkN-2(t) H — RkN{t), 

1-p 1 

RkN{t + 1) = — - — RkN-l(t) + ~RkN+l{t), (39) 

for k = 1, 2, • • • . The dynamics is different near the origin, 

i2o(t + l) = + 

Rtit + l) = R (t) + ^R 2 (t), (40) 

where 5" = YlT=i[RkN-i(t) + R^it))- At large times the system reaches a stationary state 
with RkN+i(t + 1) = RkN+i(t), and Eqs. (|58j) . (j3Uj) and (jjO)) simplify to yield for the first 
period (k — 0), 

Ri = 2R - ipS' for % = 1, 2, ■ ■ ■ , N - 1; 

= ^(2^ -A^'), (41) 

while for all other periods (k > 1) one can obtain 

RkN+i = (i + ±)RkN - i{l - p)RkN-u for i = 0, 1, JV - 1; 
R kN+N = YZ_~p [(^ + - A(l - p)R kN ^\ . (42) 

In order to solve Eqs. (|4*T|) and (J4*2*|) it is natural to look for a solution in the form 

R kN+l = 2R e ak -%B{k), (43) 
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for all k and % except the origin (k = i = 0). Note that the factor of 2 in front of Rq in 
Eq. (pfHj) is necessary to ensure the consistency between the solutions for k = from Eq. 
({H} and for k > 1 [see Eq. (|42j)]. The solutions can be obtained using the same approach 
as was discussed in the previous section for the continuous-time forward BBM. The final 
expressions are 

1 — x 



Ro 



N+ [N-p(N- l)]x 5 



RkN+i — 2i? X 



l-^(l-(l-p)ar) 



(44) 



where we again define x = e a , and it can be shown that 



x=^(a-Vc?=i), a = ^-±j-(l-p)(N-l). (45) 

The dynamic properties of the particle in this model can be calculated as presented 
above for the continuous-time forward BBM. For the velocity we obtain the following formal 
expression, 

oo 

v = Ro + lY. RkN - ( 46 ) 

k=l 

The validity of this equation can be seen from the fact that the forward and backward 
hopping rates cancel each other on all sites except at the origin and the sites labeled kN 
(k = 1,2, •••). From Eq. (}4*4"j) it follows that R^n = 2R x k . Substituting this result into 
yields 

„ = fi 1 + ' p - 1)x , (47) 
I — X 



(48) 



that can be transformed with the help of Eq. ()44|) into 

1 — (1 — p)x 



N+ [N-p(N-l)]x 
Finally, replacing x according to (|4"Hj) and using = 1/c we obtain 



where 




z = [2c +P( 2 - c)] { -i + t /i - [2c+ ^_m- c)]. > ■ < 50) 

After some rearrangement, it can be shown that this formula is identical to the expression 
for the velocity obtained by Antal and Krapivsky [7J. But note again that our derivation is 
much simpler. 
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Random distribution of bridges. 

So far we considered BBM with weak links positioned at equal distances from each other, 
i.e., periodic distribution of bridges. Our approach also allows to easily calculate the prop- 
erties of the systems with random distribution of bridges. In this case, each link can be a 
bridge with the probability c every time the particle visits it, and the link can also burn 
with the probability p. This suggests that BBM with random distribution of bridges is 
equivalent to BBM with periodic distribution for N = 1 and with the probability of burning 
equal to pc. Then all dynamic properties can be calculated as explained above by modifying 
already derived results with c — ► 1 and p — > pc. For example, the velocity of the particle in 
the continuous-time forward BBM with random distribution of bridges is obtained from Eq. 
(123) as 

V r and(c,p) = V per (l, Cp) = y/pc. (51) 

Also for the random distribution of bridges in the continuous-time forward-backward model 
the calculations based on the method outlined above produce the following expression, 

V ran d(c,p) = pc(2 - pc) . (52) 

In the limit of c 1 and p <C 1 it reduces to \J2pc, which is in exact agreement with 

n 

the result obtained for this case in the Ref. p. It should be noted again, that there are 
different possibilities of introducing randomness in the distribution of bridges that differ from 
our presentation, and that could effect the dynamic properties of the system. However, our 
theoretical approach can be applied for these cases too. 



DISCUSSIONS 



To check the validity of our theoretical approach we performed extensive Monte Carlo 
computer simulations for several BBMs. A general algorithm for continuous-time and 
discrete-time random walks, discussed in detail in Ref. 11 1, has been utilized. Veloci- 
ties have been calculated from the single simulations, consisting of 10 6 -10 7 steps, via the 
basic formula V — (x) / 1. Diffusion constants have been determined using the expression 
D = ((x 2 ) — (x) 2 )/2t, and typically more than 10 5 simulation runs have been averaged over 
to reduce the stochastic noise in the data. The results of simulations, as shown in Figs. 2, 
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4, 5A and 6, are in all cases in excellent agreement with available analytical calculations, 
supporting our analytical predictions. 

Theoretical calculations show that the velocity of the particle in BBM is a monotonously 
increasing function of the concentration of weak links and burning probability: see Figs. 2, 
4 and 5. However, the behavior of the dispersion is more complex, as can be judged from 
Figs. 2 and 6. For fixed values of p it decreases as a function of c, although note large 
fluctuations when c — > 0. Increase in the burning probability first lowers significantly the 
dispersion, and then it starts to increase slowly. Our computer simulations indicate (see Fig. 
6a) that when every link in the system is a potential bridge (c = 1) the dispersion is always 
equal to 1/2 for any p ^ 0. 

Dynamic properties of the particle in BBM depend on the distribution of bridges, as 
illustrated in Fig. 5A. Contrary to naive expectations, at the same concentration c and 
the burning probability p, the random walker moves faster for the random and not for the 
periodic distribution of bridges, although the difference decreases for lower values of p. This 
surprising observation can be explained using the following arguments. When burning of 
the bridge can take place with a significant probability, the rate-limiting step in the forward 
motion of the particle is to find the weak link. For the case of the random distribution, when 
any link can turn out to be a bridge with the probability c, the time to find the bridge is of 
the order of 1/c. In the periodic distribution of bridges the particle has to diffuse a distance 
N = 1/c before reaching the bridge, and the diffusion time is of the order of N 2 = 1/c 2 . 
Thus, the particle moves faster in the random distribution of bridges. 

Our theoretical calculations also show that details of burning process (forward versus 
forward-backward) are affecting the dynamics of the particle. Under the same conditions, 
the random walker moves faster in the forward-backward BBM, and this effect is significant 
for large p and c: see Fig. 5B. These results are expected since the burning of the bridge fuels 
the forward motion of the particle. The possibility of burning in the forward and backward 
motion accelerates the motion of the random walker. 

Developed theoretical method allows us to analyze experimental data on the motion of 
collagenases motor proteins |4j. Using fluorescence correlation spectroscopy it was deter- 
mined that the diffusion coefficient is D = 8 ±1.5 x 10~ 13 m 2 s _1 and the mean drift velocity 
is V = 4.5 ± 0.36 x 10 -6 m s _1 Q. Our theoretical results are obtained assuming that 
the lattice spacing and transition rates are equal to unity. To restore the correct units we 
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assume that the motor protein moves in discrete steps of size a. Although the exact value 
of a is unknown we can reasonably take a ~ 10 nm, which is of the same order as step sizes 
of other molecular motors Q|. From the structure of the collagen it is known that _c 
sites, that correspond to the bridges, are separated by the distance A = 300 nm 
suggests that BBM with periodic distribution of bridges and c— a/A = 1/30 is appropriate 
to describe the dynamics of the system. Since at low concentrations the predictions for the 
dynamic properties in the forward and forward-backward BBM are very similar (see Fig. 
5B), for the velocity we use Eq. (|23j) multiplied by av , 

V{c,p) = av 2p ° , (53) 

\/p {1 — cy + Ape 

where vo is the transition rate. For the diffusion coefficient analytical results are not known, 
but for cCl the limiting result of c = can be used, yielding 

D = a 2 v . (54) 

Eqs. (p)3*j) and (p>4*|) with 2 unknowns, p and Vq, can be easily solved, and we obtain 

p ~ 0.283, v ~ 8 x 10 3 s" 1 . (55) 

Our estimate of the burning probability p ~ 28.3% is larger than the result p = 10% 
obtained in Monte Carlo simulations 4]. However, the difference can be attributed to 
different versions of BBM used in our theoretical analysis and in computer simulations jj| , 
and to the uncertainty in the determination of the step size a. 

It is also interesting to analyze the efficiency of the motor protein driven by burning of 
bridges. Our theoretical analysis indicates that the motion of the particle in the system 
with burnt bridges can be viewed as a biased random walk. The dynamic properties of the 
random walker effectively can be written as 

V = a(u eff - w eff ), D = ^a 2 {u eff + w eff ), (56) 



where u e ff and w e ff are effective transition rates of the 



exerted by the molecular motor can be found from 



3,Q 



Diased random walk. Then the force 



F= kBT ln ML. (57 ) 
a w eff 
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Substituting experimental values for V and D into Eqs. (|3B)l to obtain the effective rates 
u e ff and w e ff, the exerted force can be calculated, yielding F ~ 0.023 pN, which is two 
orders of magnitude weaker than the forces produced by other motor proteins Q] • However, 
this force is exerted over the distance of A = 300 nm, which is 10-100 times larger than 
for conventional molecular motors, producing W = F x A ~ 6.9 pNxnm of work. The 
cleavage of one peptide bond, that fuels the motion of collagenase, liberates E ~ 20.9 
kJ/mole or 34.8 pNxnm of energy [4]. The calculated efficiency of the collagenase motor 
protein, v = W/E ~ 19.8%, is only slightly larger than the efficiency of 15% obtained from 
the experimental estimates [4J. Thus the collagenase is a weak molecular motor and its 
efficiency is not very high in comparison with other motor proteins. 



SUMMARY AND CONCLUSIONS 



A general theoretical analysis of particle dynamics in different one-dimensional burnt- 
bridge models is presented. It is shown that when the probability of burning is equal 
to unity, the system can be mapped into the model of the single particle hopping along 
the periodic infinite lattice. This allows us to obtain all dynamic properties in exact and 
explicit form. For general case of p < 1 we developed a method of analyzing the dynamics 
by considering a reduced kinetic scheme with respect to the last burnt bridge. It is found 
that distribution of weak links and the details of the burning mechanism strongly influence 
the particle dynamics in BBM. Theoretical calculations are fully supported by Monte Carlo 
computer simulations. Application of this theoretical method to experiments on collagenase 
motor proteins provides the estimates of the burning probability and the exerted force, 
leading to conclusion that this molecular motor is rather weak and not very efficient, in 
contrast to conventional motor proteins. 

Although the presented theoretical approach investigates explicitly the dynamic proper- 
ties of many BBMs and it is successfully applied for analysis of experimental data, there are 
several problems that should be addressed in the future studies. We have not been able to 
obtain expressions for the diffusion coefficient in the general case of p < 1. It is also not 
clear what is the origin of jump in the dispersion in the limit of c — > 0. It will be interesting 
to investigate the motion of the motor protein on periodic lattices [h| with bridge burning, 
and it is reasonable to suggest that the approach developed in this paper can be generalized 
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for this case too. In addition, a more realistic approach will be to allow the molecular motor 

n 

to jump between parallel tracks |4J. Future experiments on single motor proteins that utilize 
bridge burning mechanisms will be important to test the validity of the presented theoretical 
approach. 



ACKNOWLEDGMENTS 

The authors would like to acknowledge the support from the Welch Foundation (grant C- 
1559), the U.S. National Science Foundation (grant CHE-0237105) and Hammill Innovation 
Award. 



[1] H. Lodish et. al., Molecular Cell Biology, (W.H. Freeman and Company, New York, 2000). 
[2] J. Howard, Mechanics of Motor Proteins and Cytoskeleton, (Sinauer Associates, Sunderland 
Massachusetts, 2001). 

[3] D. Bray, Cell Movements. From Molecules to Motility, (Garland Publishing, New York, 2001). 
[4] S. Saffarian, I.E. Coller, B.L. Marnier, E.L. Elson, and G. Goldberg, Science 306, 108 (2004). 
[5] S. Saffarian, H. Qian, I.E. Coller, E.L. Elson, and G. Goldberg, Phys. Rev. E 73, 041909 
(2006). 

[6] J. Mai, I.M. Sokolov, and A. Blumen, Phys. Rev. E 64, 011102 (2001). 

[7] T. Antal and PL. Krapivsky, Phys. Rev. E 72, 046104 (2005). 

[8] A.B. Kolomeisky and M.E. Fisher, J. Chem. Phys. 113, 10867 (2000). 

[9] B. Derrida, J. Stat. Phys. 31, 433 (1983). 
[10] A.B. Kolomeisky and M.E. Fisher, to appear in Ann. Rev. Phys. Chem, (2007). 
[11] J. F. McCarthy, J. Phys. A: Math. Gen. 26, 2495 (1993). 

[12] M.E.Fisher and A.B. Kolomeisky, Proc. Natl. Acad. Sci. USA 96, 6597 (1999). 



17 



Figure Captions: 

Fig. 1. A schematic view of the motion of the particle in the system with burnt bridges. 
Thick solid lines are strong links, while thin solid lines are weak links, or "bridges" that can 
be burnt. Dotted lines are already burnt bridges. The particle can jump with equal rates 
to the right or to the left unless the link is broken. 

Fig. 2. Dynamic properties of continuous-time burnt-bridge model with periodic distribution 
of bridges for p = 1. Lines are from analytic calculations, while symbols are obtained from 
Monte Carlo computer simulations - see text for details. 

Fig. 3. Reduced kinetic schemes for continuous-time systems with periodic distribution of 
bridges: a) forward BBM (only transition rates not equal to unity are shown); b) forward- 
backward BBM (only transition rates not equal to 1/2 are shown). The origin corresponds 
to the right end of the last burnt bridge. 

Fig. 4. A) Velocity as a function of concentration of bridges for different burning probabil- 
ities for the periodic forward continuous-time BBM. Solid lines are results of analytical cal- 
culations, and symbols are from Monte Carlo simulations. Curve a) corresponds to p — 0.9; 
curve b) corresponds to p — 0.7; curve c) corresponds to p — 0.5; and curve d) corresponds 
top = 0.1. 

B) Velocity as a function of burning probability for different concentrations of weak links 
for the periodic forward continuous-time BBM. Solid lines are results of analytical calcula- 
tions, and symbols are from Monte Carlo simulations. Curve a) corresponds to c = 1; curve 
b) corresponds to c = 1/2; curve c) corresponds to c = 1/3; and curve d) corresponds to 
c= 1/10. 

Fig. 5. A) Velocity as a function of concentration for different distributions of bridges for the 
periodic forward continuous-time BBM. Solid lines are results of analytical calculations, and 
symbols are from Monte Carlo simulations. Curve a) corresponds to random distribution 
with p — 1; curve b) corresponds to periodic distribution with p — 1; curve c) corresponds 
to random distribution with p = 0.5; and curve d) corresponds to periodic distribution with 
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p = 0.5. 

B) Velocity as a function of concentration for different burning realizations for the periodic 
continuous-time periodic BBM. Solid lines are results of analytical calculations. Curve a) 
corresponds to forward-backward BBM with p = 0.8; curve b) corresponds to forward BBM 
with p = 0.8; curve c) corresponds to forward-backward BBM with p = 0.2; and curve d) 
corresponds to forward BBM with p = 0.2. 

Fig. 6. Dispersion in the periodic forward continuous-time periodic BBM from Monte Carlo 
simulations: A) as a function of burning probability; and B) as a function of concentration 
of bridges. 
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Fig.3b 
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Fig.4A 



24 




Fig.4B 
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Fig.5A 
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Fig.5B 
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Fig.6A 
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Fig.6B 



